Preparation and detection of magnetic quantum phases in optical superlattices 
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We describe a novel approach to prepare, detect and characterize magnetic quantum phases in ultra-cold 
spinor atoms loaded in optical superlattices. Our technique makes use of singlet-triplet spin manipulations in an 
array of isolated double well potentials in analogy to recently demonstrated quantum control in semiconductor 
quantum dots. We also discuss the many-body singlet-triplet spin dynamics arising from coherent coupling 
between nearest neighbor double wells and derive an effective description for such system. We use it to study 
the generation of complex magnetic states by adiabatic and non-equilibrium dynamics. 
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Recent advances in the manipulations of ultra-cold atoms 
in optical lattices have opened new possibilities for exploring 
complex many -body systems jl]]. A particular topic of con- 
tinuous interest is the study of quantum magnetism in spin 
systems JUSS]- By loading spinor atoms in optical lattices it 
is now possible to "simulate" exotic spin models in controlled 
environments and to explore novel spin orders and phases. 

In this Letter we describe a new approach for prepara- 
tion and probing of many-body magnetic quantum states that 
makes use of coherent manipulation of singlet-triplet pairs of 
ultra-cold atoms loaded in deep period-two optical superlat- 
tices. Our approach makes use of a spin dependent energy 
offset between the double-well minima to completely control 
and measure the spin state of two-atom pairs, in a way anal- 
ogous to the recently demonstrated manipulations of coupled 
electrons in semiconductor double-dots 15[]. As an example, 
we show how this technique allows one to detect and analyze 
anti-ferromagnetic spin states in optical lattices. We further 
study the many-body dynamics that emerge when tunneling 
between nearest neighbor double wells is allowed. As two 
specific examples, we show how a set of singlet atomic states 
can be evolved into singlet-triplet cluster-type states and into a 
maximally entangled superposition of two anti-ferromagnetic 
states. Finally, we discuss the use of our projection technique 
to probe the density of spin defects (kinks) in magnetic states 
prepared via equilibrium and non-equilibrium dynamics. 

The key idea of this work is illustrated by considering a pair 
of ultra-cold atoms with two relevant internal states, which 
we identify with spin up and down a =f , J, in an isolated 
double well (DW) potential as shown in Fig.l. By dynam- 
ically changing the optical lattice parameters, it is possible 
to completely control this system and measure it in an arbi- 
trary two-spin basis. For concreteness, we first focus on the 
fermionic case. The physics of this system is governed by 
three sets of energy scales: i) the on-site interaction energy 
U = Ufi between the atoms, ii) the tunneling energy of the 
cr species: J CT , and iii) the energy difference between the two 
DW minima, 2A CT for each of the two species. The a index 
in J and A is due to the fact that the lattice that the \ and 
j atoms feel can be engineered to be different by choosing 



laser beams of appropriate polarizations, frequencies, phases 
and intensities. In the following we assume that the atoms are 
strongly interacting, U J a , and that effective vibrational 
energy of each well, hu>o, is the largest energy scale in the 
system Hluq ^> U, A CT , J a , i.e deep wells. 

Singlet \s) and triplet \t) states form the natural basis for 
the two-atom system. The relative energies of these states 
can be manipulated by controlling the energy bias A a be- 
tween the two wells. In the unbiased case (U ^S> 2A CT ) only 
states with one atom per site (1, 1) are populated, as the large 
atomic repulsion energetically suppresses double occupancy 
(here, labels (m, n) indicate the integer number of atoms in 
the left and right sites of the DW). For weak tunneling and 
spin independent lattices (Jj = Jy = J, Af = A; = A) 
the states (1, l)|s) and (1, are nearly degenerated. The 
small energy splitting between them is ~ 4J 2 /U, with the 
singlet being the low energy state (Fig. la). As A is increased 
the relative energy of doubly occupied states (0, 2) decreases. 
Therefore, states (1, l)|s) and (0, 2)|s) will hybridize. When 
2A > U the atomic repulsion is overwhelmed and conse- 
quently the (0,2)|s) becomes the ground state. At the same 
time, Pauli exclusion results in a large energy splitting hojo 
between doubly occupied singlet and triplet states as the lat- 
ter must have an antisymmetric orbital wave function. Hence, 
(1, l)|t) does not hybridize with its doubly occupied counter- 
part, and its relative energy becomes large as compared to the 
singlet state. Thus the energy difference between singlet and 
triplet states can be controlled using A. 

Further control is provided by changing J a and A CT in spin 
dependent lattices (see Fig. lb). Specifically, let us now con- 
sider the regime 2A CT <C U in which only (1, 1) subspace is 
populated. Within this manifold we define ||6(] 

l*> = ^10)^-^(114) -| It)), (i) 
= tt|o> = 4=(IU> + UT», (2) 
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\t x ) = 4|o) ee ^=(| TT) - | U)), (3) 
\ty) = *»^(in> + lll» (4) 

Here v a and s are operators that create triplet and sin- 
glet states from the vacuum |0) (state with no atoms). 
They satisfy bosonic commutation relations and the constrain 
(^ tt=I y z fi a t a ) + s^s = 1, due to the physical restriction 
that the state in a double well is either a singlet or a triplet. In 
the rest of the letter we will omit the label (1, 1) for the singly 
occupied states. 

When A CT depends on spin, i.e T = A-f — Aj, ^ 0, the 
\t z ) component mixes with |s)(see Fig.lc). Note that on the 
other hand \t XlV ) remain decoupled from \t z ) and \s) . As a 
result the states \s) and \ t z ) form an effective two-level system 
whose dynamics is driven by the Hamiltonian: 

H( = -C(s f s - i\t z ) - TS Z + const, (5) 

Here £ = 2J^JJU, is the exchange coupling energy (with 
U ee u2 -( A T j + A i) 2 - ) and s z = sH z + t\s. If T = 0, ex- 
change dominates and \s) and \t z ) becomes the ground and 
first excited states respectively. However if T 3> f, exchange 
can be neglected and the ground state becomes either | ||) or 
||t) depending on the sign of T. 




FIG. 1 : (color online) a) Energy levels of fermionic atoms in a spin 
independent double well as A/U is varied: While in the regime 
2A < U, (1, l)|s) is the lowest energy state, when 2A > U, 
(0, 2)|s) becomes the state with lowest energy, b) In spin depen- 
dent potentials the two species feel different lattice parameters c) 
Restricted to the (1,1) subspace T acts as an effective magnetic field 
gradient and couples the |s) and \ t z ) states . 

These considerations indicate that it is possible to perform 
arbitrary coherent manipulations and robust measurement of 



atom pair spin states. The former can be accomplished by 
combining time-dependant control over £, T to obtain effec- 
tive rotations on the spin- 1/2 Bloch sphere within \s) — \t z ) 
state. In the parameter regime of interest, £, T, can be var- 
ied independently in experiments. In addition, by applying 
pulsed (uniform) magnetic fields it is possible to rotate the 
basis, thereby changing the relative population of the \t x ^ V:Z ) 
states. Atom pair spin states can be probed by adiabatically 
increasing A until it becomes larger than {7/2, in which case 
atoms in the \s) will adiabatically follow to (0, 2)|s) while the 
atoms in \t a ) will remain in (1,1) state (Fig. la). A subse- 
quent measurement of the number of doubly occupied wells 
will reveal the number of singlets in the initial state. Such 
a measurement can be achieved by efficiently converting the 
doubly occupied wells into molecules via photoassociation or 
using other techniques such as microwave spectroscopy and 
spin changing collisions J3]. Alternatively, one can continue 
adiabatically tilting the DW until it merges to one well. In 
such a way the \s) will be projected to the (0, 2)|s), while the 
triplets will map to (0, 2)\t a ). As (0, 2)\t a ) has one of the 
atoms in the first vibrational state of the well, by measuring 
the population in excited bands one can detect the number of 
initial \t a ) states. Hence the spin-triplet blockade |5f| allows 
to effectively control and measure atom pairs. 

Detection and diagnostics of many-body spin phases such 
as antiferromagnetic (AF) states is an example of direct ap- 
plication of the singlet-triplet manipulation and measurement 
technique. The procedure to measure the AF state population 
is the following; after inhibiting tunneling between the various 
DWs, one can abruptly increase T, such that the initial state 
is projected into the new eigenstates | ||) and | |f) at time 
t = tq. For r > tq T can then be adiabatically decreased to 
zero, in which case the | f|) pairs will be adiabatically con- 
verted into \s) and | ||) pairs to \t z ). Finally, the singlet pop- 
ulation can be measured using the spin blockade. As a result, 
a measure of the doubly occupied sites (or excited bands pop- 
ulation) will detect the number of | ||) pairs and thus probe 
antiferromagnetic states of the type | till ••■)■ 

These ideas can be directly generalized to perform mea- 
surements of the more complex magnetic states that can be 
represented as products of two atom pairs. For example, a 
pulse of RF magnetic field can be used to orient all spins, 
thus providing the ability to detect | AF) states aligned along 
an arbitrary direction. Moreover, one can determine the rela- 
tive phase between singlet and triplet pairs in | AF) states of 
the form Yl \ s ) + e \tz) by performing Ramsey-type spec- 
troscopy. After letting the system evolve freely (with T = 0) 
so that the \s) and \t z ) components accumulate an additional 
relative phase due to exchange, a read-out pulse (controlled 
by pulsing T) will map the accumulated phase onto popula- 
tion of singlet and triplet pairs. To know (f> is important as it 
determines the direction of the anti-ferromagnetic order. Fur- 
thermore, by combining the blockade with noise correlation 
measurements JH it is possible to obtain further information 
about the magnetic phases. While the blockade probes local 
correlation in the DWs, noise measurements probe non-local 
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spin-spin correlations and thus can reveal long range order. 

Before proceeding we note that similar ideas to that out- 
lined above can be used forbosonic atoms if initially no \t x>y ) 
states are populated. The latter can be done by detuning the 
\t x>y ) states by means of an external magnetic field. In the 
bosonic case the doubly occupied t z states will be the ones 
that have the lowest energy. They will be separated by an 
energy Hluq from the doubly occupied singlets as the latter 
are the ones that have antisymmetric orbital wave function in 
bosons. Consequently, the role of \s) in fermions will be re- 
placed by \t z ) in bosons. The read-out procedure would then 
be identical to that described above, while the coherent dy- 
namics will be given by the Hamiltonian Eq.© apart from 
the sign change £ — > — £. 

Up to now our analysis has ignored tunneling between dif- 
ferent DWs, but in practice this tunneling can be controlled by 
tuning the lattice potential. How will singlet and triplet pairs 
evolve due to this coupling? We will now discuss the many- 
body dynamics that emerges when nearest neighbor DW tun- 
neling is allowed, i.e. t a > 0. When atoms can hop between 
DWs, the behavior of the system will depend on the dimen- 
sionality. For simplicity we will restrict our analysis to a ID 
array of N double-wells, where t a corresponds to hopping en- 
ergy of ex-type atoms between the right site of the j th — DW 
and the left site of the (j + l)*' 1 - DW. 

In the regime Ja,t a , A CT <C U, multiply occupied wells 
are energetically suppressed and the effective Hamiltonian is 
given by H e ** ~ Hj+H t ■ Here the first term corresponds to 
the sum over N independent Hj Hamihonians (see Eq.©), 

Hj = Y^,j=i Hj> eacrl °f which acts on its respective j th - 
DW. On the other hand H t is non-local as it couples different 
DWs and quartic as it consists of terms with four singlet-triplet 
operators |9[]. The coupled DWs system is in general complex 
and the quantum spin dynamics can be studied only numeri- 
cally. However, there are specific parameter regimes where an 
exact solution can be found. For this discussion we will set 
A ff = 0. If if/i; — > 0, and at time r = 0, no \t x ), \t y ) triplet 
states are populated, their population will remain always zero. 
Consequently, in this limit, the relevant Hilbert space reduces 
to that of an effective spin one-half system with |s) and \t z ) 
representing the effective ±1/2 states, which we denote as 
| ft) and | 4J.). H t couples such effective spin states. In the 
restricted Hilbert space H e ^ maps exactly to an Ising chain 
in a magnetic field: 

3 3 

where a a are the usual Pauli matrices which act of the effec- 
tive | ft) and | JJ) spins. In terms of singlet-triplet operators 
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they are given by d] = (s}sj - V ZJ t ZJ ), a* - Oj o xj -r v zj0j 

t - -+ t 2 t 2 

andaj = {Qt g j — V X jSj)/i. Here A z = 5^7 - xjjj and the up- 
per and lower signs are for fermions and bosons respectively. 
For fermions in the lowest vibrational level the onsite interac- 
tion energy between the same type of atoms U^,Un — > oo 
due to the Pauli exclusion principle. 



The ID quantum Ising model exhibits a second order quan- 
tum phase transition at the critical value \g\ = |A 2 /C| = 1. 
For fermions (upper sign) when g <C 1 the ground state cor- 
responds to all effective spins pointing up, i.e \G) = | ft . . . ft 
) = Hj|s)j. On the other hand when g ^> 1, there are two 
degenerate ground states which are, in the effective spin ba- 
sis, macroscopic superpositions of oppositely polarized states 
along x. In terms of the original fermionic spin states this su- 
perposition correspond to the states lAi*^) = ^=(| Tl • ■ • T-L 
) ± | |t • ■ • IT))- Therefore, by adiabatic passage one could 
start with \G) and convert it into AF state(s). Due to vanish- 
ing energy gap at the quantum critical point g = L adiabatic- 
ity is difficult to maintain asiV^TOmElElB- In 
that respect, our projection scheme is useful to test adiabatic 
following. It can be done either by measuring the number of 
| 11} pairs in the final state or by adiabatically ramping down 
g back to zero and measuring the number of singlets/triplet 
pairs. The remaining number of triplets will determine the 
number of excitations created in the process. 

We now turn to non-adiabatic dynamics. We will discuss 
the situation where initially the system is prepared in a prod- 
uct of singlet states (A 2 = ground state ) and then one lets it 
evolve for r > with a fixed |A Z | > 0. Generically the cou- 
pling between DWs results in oscillations between singlet and 
triplet pairs with additional decay on a slower time scale. We 
present two important special cases involving such dynamics: 

i) Singlet-triplet cluster state generation: If the value of X z 
is set to be |A 2 | £, then the Hamiltonian reduces to a pure 
Ising Hamiltonian and thus at particular times, r c , given by 
A 2 T c /ft = 7r/4 mod 7r/2 the evolving state becomes a d — 1 
cluster state \C) in the effective spin basis 11511 . Up to single 
spin rotations \C) = ^®* =1 {\ ft)^| +1 + | ^) 3 )- Cluster 
states are of interest for the realization of one-way quantum 
computation proposals where starting from the state \C) com- 
putation can be done via measurements only. Preparation of 
cluster states encoded in the logical ft, JJ. qubits may have sig- 
nificant practical advantages since the ft, JJ. states have zero 
net spin along the quantization axis and hence are not affected 
by global magnetic field fluctuations. Additionally, the use 
of such singlet-triplet states for encoding might allow for the 
generation of decoherence free subspaces insensitive to col- 
lective and local errors lfl6ll and for alternative schemes for 
measured-based quantum computation lfl7ll . 

ii) Non-equilibrium generation and probing of AF corre- 
lations: The second situation is when the value of X z is set 
to the critical value, |A 2 | = £ (or g = 1). We will first 
focus on the fermionic system A 2 > 0. To discuss it, we 
remind that the dynamics driven by is exactly solv- 
able as H e ff can be mapped via the Jordan Wigner trans- 
formation into a quadratic Hamiltonian of fermionic opera- 
tors which can be diagonalized by a canonical transforma- 
tion ifiol Using such transformation it is possible to 
show that at specific times, the shortest of them we denote 
by r m « H^^, long range AF correlations build up and for 
small atom number the state approaches \AF + ). To quan- 
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tify the resulting state in Fig. 2(inset) we plot the fidelity, 
defined as .Fi(t to ) = |(AF + |^(r m )) ff= i| 2 , as a function of 
N. The figure shows that while an almost perfect | AF+) is 
dynamically generated for small TV, its fidelity exponentially 
degrades with increasing atom number. 

However, the fidelity is a very strict probe, as it drops to 
zero when a single spin is flipped. As TV increases the sys- 
tem ends at r m in a quantum superposition of states like 
| • • • =><S=<=<=-^=-^==>=>=>=>=>=><^= • ■ ■ ) with finite do- 
mains of "effective spins" pointing along ±x, separated by 
kinks where the polarization of the spins change its orien- 
tation (we used the convention | fj.) = | =>}). Conse- 
quently, one gets more realistic information about the AF 
order of the state, by measuring the average size of the do- 
mains or the average density of kinks, the latter defined as 



1 

2N 



Ei(l-Wr)|5f5f +1 |V(T))) 



Our read-out technique can be used to detect the kink- 
density as for an arbitrary fixed g energy conservation imposes 
a relation between v and the triple t-z density, N t : 



1 N t (T,g) 

= ~ ■ 

2 g 



(7) 



A simple analytical expression for Nt(r,g) can be ob- 
tained by using the Jordan Wigner transformation fioll : 



N t (T,g) 



±\ 2 



-1 sin 2 (27rfc/jV) sin 2 (2u fc -r) 




where fiwi 



C\/g 2 + 1 + 2gcos(2irk/N) are quasi-particle frequencies of 
H e ff. The fact that it remains always below 0.2 (see Fig. 
2) confirms the idea that regardless of the reduced fidelity at 
large N, the state does retain AF correlations. We point out 
that \AF + ) states are only generated at g = 1, a feature that 
illustrates the special character of the critical dynamics. 
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FIG. 2: Using the the Jordan- Wigner transformation 1 1 Oil we cal- 
culated the density of kinks vs N at r = r m and the fidelity 
| (^(t to ) I AF + ) | 2 vs N (inset). Our projection technique can be used 
to measure v(t) as it is directly related to the triplet density, Nt(r) 
(see Eq.(7)). 

Let us now discuss the bosonic case. If X z > 0, the 
fermionic results apply for bosons by simply interchanging 
the role of \s) <-» \t z ). On the other hand if X z < 0, not only 
one has to interchange |s) «-> \t z ) but additionally, the adia- 
batic and non-equilibrium dynamics will generate, instead of 



macroscopic superpositions of AF states along the x-direction 
in the effective spin basis. With these modifications, the re- 
sults derived for fermions hold for bosons 12 ill . 

Before concluding we briefly mention that spin dependent 
superlattices of the form 

V = ( A J + Bj<7 z )cos 2 [kz/j + 9j] (8) 

.7=1,2 

can be experimentally realized by superimposing two inde- 
pendent lattices, generated by elliptic ally polarized light, one 



with twice the periodicity of the other Ill8l 20ll . Com- 
plete control over the DW parameters is achieved by control- 
ling the phases (which determine A), intensities (which deter- 
mine U, J and t) and polarization of the laser beams (which 
allow for spin dependent control). For example lattice con- 
figurations with tf <C t± can be achieved by setting the laser 
parameters such that B\ = and A<z = B2 3> 1. 

In summary we have described a technique to prepare, de- 
tect and manipulate spin configurations in ultra-cold atomic 
systems loaded in spin dependent period-two superlattices. 
By studying the many-body dynamics that arises when tun- 
neling between DWs is allowed, we discussed how to dynam- 
ically generate singlet-triplet cluster states and AF cat states, 
which are of interest for quantum information science, and 
how to probe AF correlations in far from equilibrium dynam- 
ics. Even though in this Letter we restrict our analysis to ID 
systems the ideas developed here can be extended to higher 
dimensions and more general kinds of interactions. 
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In this case a different sign in the definition of kink density 
of of+i — > — &f&f +1 is required. 



